arXiv: 1505.02189vl [cond-mat.quant-gas] 8 May 2015 


Enhanced stripe phases in spin-orbit-coupled Bose-Einstein condensates in ring 

cavities 


Farokh Mivehvar and David L. Feder* 

Institute for Quantum Science and Technology, Department of Physics and Astronomy, 
University of Calgary, Calgary, Alberta, Canada T2N INf 
(Dated: May 12, 2015) 

The coupled dynamics of the atom and photon fields in optical ring cavities with two counter- 
propagating modes give rise to both spin-orbit interactions as well as long-ranged interactions be¬ 
tween atoms of a many-body system. At zero temperature, the interplay between the two-body and 
cavity-mediated interactions determines the ground state of a Bose-Einstein condensate. In this 
work, we find that cavity quantum electrodynamics in the weak-coupling regime favors a stripe-phase 
state over a plane-wave phase as the strength of cavity-mediated interactions increases. Indeed, the 
stripe phase is energetically stabilized even for condensates with attractive intra- and inter-species 
interactions for sufficiently large cavity interactions. The elementary excitation spectra in both 
phases correspond to linear dispersion relation at long wavelengths, indicating that both phases 
exhibit superfluidity, though the plane-wave phase also displays a characteristic roton-type feature. 
The results suggest that even in the weak coupling regime cavities can yield interesting new physics 
in ultracold quantum gases. 


I. INTRODUCTION 

The experimental realization of Bose-Einstein conden¬ 
sation (BEC) has opened many opportunities for real¬ 
izing new many-body phases [1-3]. Ultracold atoms 
trapped in laser-generated optical lattice potentials ex¬ 
perience crystalline environments and exhibit a variety of 
intriguing phenomena [4], most notably the superfluid- 
Mott-insulator phase transition [5]. There are numerous 
proposals for inducing gauge fields in quantum gases by 
means of laser light [6], and recently abelian [7] and non- 
abelian [8] gauge fields have been realized. In the latter 
work an equal combination of Rashba and Dresselhaus 
spin-orbit (SO) couplings were induced via two-photon 
Raman transitions. These developments have set the 
stage for realizing topological states in these systems [9]. 

The single-particle energy dispersion of a SO-coupled 
atom is a momentum-space double well, which is two-fold 
degenerate in the symmetric case [7]. In a Bose-Einstein 
condensate (BEC) of atoms, the two-body interactions 
lift this degeneracy and drive the BEC into either a plane 
wave phase (PWP) or a stripe phase (SP), depending on 
the strength and sign of the intra- and inter-species two- 
body interactions [10-13]. In the PWP, all atoms con¬ 
dense into one of the two single-particle energy minima, 
while the SP is a superposition state of the minima and 
the total BEC density exhibits faint fringes [14]. Addi¬ 
tional phases are found for fully three-dimensional SO 
interactions [15]. When a SO-coupled quantum gas is 
confined in an optical lattice, the ground state of the 
system exhibits a variety of magnetic orderings in the 
Mott-insulator regime, such as ferromagnetic, antiferro¬ 
magnetic, spin spiral, vortex and antivortex crystals, and 
skyrmion crystal phases [16-18]. The superfluid to Mott- 
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insulator phase transition of SO-coupled quantum gases 
has also been investigated [16, 19]. 

In laser-based approaches to generating SO couplings, 
the radiation field is treated classically and one ignores 
the back-action of the atoms on it. Confining the radia¬ 
tion field to within an optical cavity leads to a coherent 
exchange of energy and momentum between atoms and 
photons [20]. The back-action of the atoms on the photon 
fields is no longer negligible, leading to complex coupled 
dynamics of the matter and radiation fields in which both 
entities are affected by one another and must be treated 
on the same footing [21]. As a consequence, cavity- 
mediated long-range interactions are induced between 
atoms, yielding novel collective phenomena in atomic 
systems [22]. A few schemes have been recently pro¬ 
posed to induce SO coupling in ultracold atoms via cavity 
quantum electrodynamics [23-26] and to couple a laser- 
induced SO-coupled BEC to the cavity field [27]. These 
schemes exhibit a wealth of physics, including strong syn¬ 
thetic magnetic fields, a cavity-mediated Hofstadter spec¬ 
trum, and a variety of magnetic orders. 

In this work we investigate the ground state and the 
elementary excitations of a spinor BEC at zero tempera¬ 
ture subject to ring-cavity-induced SO interactions [23]. 
Here we consider lossy cavities where a steady-state pho¬ 
ton population is maintained by the application of ex¬ 
ternal pump lasers. The cavity photons mediate infinite- 
range interactions between atoms, whose strengths can 
be tuned experimentally by adjusting the amplitudes of 
the pump lasers. The sign of these interactions can be 
made positive or negative depending on the cavity detun¬ 
ing, the frequency difference between the applied pump 
lasers and the cavity. These cavity-mediated interac¬ 
tions compete with the inherent two-body interactions 
between atoms to determine the ground state of the SO- 
coupled BEC. In particular, stripe phases are always fa¬ 
vored when positive cavity-mediated interactions domi¬ 
nate the two-body-interactions, even in the case where 
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the intrinsic atomic interactions (both intra- and inter¬ 
species) are attractive. Asymmetry in the strength of 
cavity-mediated interactions for different spin compo¬ 
nents yields stripe-phase states with an arbitrary num¬ 
ber of atoms in the left or right minimum of the single¬ 
particle dispersion relation, so that the magnetization 
varies continuously from zero in the stripe phase to unity 
in the plane-wave phase. This behavior allows us to iden¬ 
tify a novel stripe-phase order parameter, and to identify 
its associated mean-field critical exponent. 

Consideration of the quantum fluctuations around the 
mean-field ground states reveals that the particle-hole el¬ 
ementary excitation spectra in both PWP and SP have 
the usual linear sound-like dispersion relation at long 
wavelengths, an indication of superfluidity. In the PWP, 
the dispersion relation also exhibits a roton-type feature 
at the same wave vector that charactizes the fringe peri¬ 
odicity in the SP, which could be used experimentally as 
a distinguishing feature. The critical transition between 
the PWP and SP occurs when the energy of this mini¬ 
mum falls below zero. Unlike for the PWP, in the SP the 
speed of sound depends strongly on the cavity-mediated 
interactions. The speed of sound is found fall below zero 
at a critical value of the cavity interactions and inter¬ 
species interactions strength, but this appears to signal 
a phase transition to a phase-separated state. Overall, 
the ring-cavity environment provides an experimentally 
convenient framework for exploring exotic ground states 
of SO-coupled BECs. 

The manuscript is organized as follows. In Section II, 
we start from the full atom-photon Hamiltonian den¬ 
sity for a lossy but pumped cavity, to derive an effective 
atomic Hamiltonian with the photon fields eliminated. 
The ground state of this effective Hamiltonian is explored 
in Section HI using both a variational method and by 
solving the generalized Gross-Pitaevskii equations. The 
remainder of this Section is devoted to an analysis of the 
elementary excitations. A discussion of the results and 
conclusions are found in Sec. IV. Appendices A and B 
provide details of the adiabatic elimination of the atomic 
excited state and cavity fields, respectively. 


II. MODEL AND HAMILTONIAN 

Consider spin-1 bosonic atoms inside a ring cavity with 
two driven counter-propagating running modes 
and where Aj is the annihilation operator for 

the photon in jth mode with wave vector kj = LOjjc 
and 2 ; is the direction along the cavity axis. Without 
loss of generality, one can assume that the wave vec¬ 
tors ki and ^2 of the two modes are approximately equal 
to each other, kR = ki ~ ^2 [28]. The mode 
(A 2 e“*^^^) propagates to the right (left) and solely in¬ 
duces the atomic transition |a) —>■ |e) (|5) —)• |e)), where 
{|a), | 6 )} are non-degenerate pseudospin states of interest 
and |e) is an excited state. The two cavity modes Aj are 
assumed to be sufficiently populated to justify omitting 


associated degenerate modes A'. In principle, a state- 
independent external potential I 4 xt(r) would need to be 
imposed to confine atoms inside the cavity. The single¬ 
particle Hamiltonian density in the dipole and rotating- 
wave approximations is 

H«=HW+ifcav+H«, (I) 


with 


^at 




^3x3 


E 


S-r(T-T 


T^{a,b,e} 

Hca-v = ^ E E _ H.c.) , 


i=i.2 


i=i.2 


A^aeb) + H.C. 


where Et are the internal atomic-state energies, arr' = 
jr) (r'l, and / 3 X 3 is the identity matrix in the internal 
atomic-state space. The atom-photon coupling for the 
transition t ^ t' is denoted and H.c. stands for 

the Hermitian conjugate. The cavity mode Aj is driven 
by a pump laser with frequency ujpj and amplitude rjj, 
indicated by the second sum in iLcav In this work, in or¬ 
der to simplify the analytical calculations, I 4 xt(r) is set 
to zero. In reality, one might imagine a very weak (al¬ 
most unbound) confining potential along the cavity axis 
z but a standard harmonic trap in the radial direction. 
The details of the transverse confining potential are not 
important for the analysis presented in this work. 

After expressing Hamiltonian (1) in the rotating frame 
of the pump lasers [29] and assuming that the atomic 
detunings Ai = uji — Sealk and A 2 = W 2 — Seb/k are 
large compared to Sba/f^ = {^b — one can adi- 

abatically eliminate the atomic excited state to obtain 
an effective Hamiltonian HgQ for the ground pseudospin 
states {[1), [2)} = {j 6 ), ja)}. The details are presented 
in Appendix A. In the limit of a very weak confining 
potential along the cavity axis z, one can assume that 
the momentum pz = hkz is a good quantum number. 
Alternatively one could consider approximately uniform 
quantum gases in a box potential where I 4 xt(r) = 0 ex¬ 
cept at the boundaries; such a potential has recently been 
realized experimentally [30]. One can then transform 
to the co-moving frame of the cavity modes by apply¬ 
ing the unitary transformation ^2 = (where 

CTz = CTii — (T 22 is the third Pauli matrix, see also Ap¬ 
pendix A). The kinetic-energy part of the Hamiltonian 
density = '^ 2 ’Hgo ^^ 2 ^ associated with the momen¬ 

tum pz, Eq. (A5), then takes the familiar form of an equal 
Rashba-Dresselhaus SO coupling: i^{pzl 2 x 2 + kknazA, 
which is characterized by a double-well energy disper¬ 
sion [ 8 ]. 

In the presence of dissipation, such as when the decay 
rate k of both cavity modes is non-zero, one should in 
principle numerically solve the associated master equa¬ 
tion [31]. That said, in the weak-coupling regime when 
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K is the dominant energy scale, k ^ ^^ae^'^he)^ the mas¬ 
ter equation approach is equivalent to including dissipa¬ 
tion in the Heisenberg equations of motion for the cavity 
fields: dtAj = —i[Aj,'H'^Q^]/h— nAj [21]. The cavity 
fields quickly reach steady states, allowing them to be 
adiabatically eliminated. Setting dtAj = 0 one obtains 
steady-state expressions for Aj that can be substituted 
into ^ 30 ^ to yield an effective atomic Hamiltonian; the 
details are relegated to Appendix B. 

The resulting effective many-body Hamiltonian reads 
+ ^51^1 + ^52^2 + 5i2hin2^ 

+ ^ UrNl + U±S+S- +U^S-S+ +2UdsNS^, 

r^l.2 

( 2 ) 



where ^(r) = (i/)i(r), ■i/) 2 (r))^ are the bosonic field opera¬ 
tors obeying the commutation relation [^/^.^(r), = 

5,-.,./(5(r — r'). Nr = J hr{'r)d^r = f 'ipl(r)'i/ir(r)d^r is the 
total atomic number operator for pseudospin r € { 1 , 2 }, 
N = JVi + JV 2 is the total atomic number operator, and 
the x-component of the total spin operator is defined 
in a usual way Sx = 5 (< 5 + + S-) using the collective 
pseudospin raising and lowering operators 3+ = = 

J V>|(r)'i/’ 2 (r)(i^r. The atoms in this system experience 
two kinds of interactions, reflected in the effective Hamil¬ 
tonian Hgg: the standard two-body contact interactions 
and the cavity-mediated long-ranged interactions. Here 
gr = Qtt denotes the two-body intra-species interaction 
strength and gi 2 the two-body inter-species interaction 
strength. The strength of the cavity-mediated interac¬ 
tions are found in Appendix B: 


r, 4h%^A,{Al-3K^) 2 

^1(2) A2(A2 + «;2)3 '?2(1): 


4h%^A, 

^±(=F) = 




A2(A2 -^^2)3 


(A2 - 3^2) 

A2(A2+k2)3 


’ 7 i( 2 )i 


(3) 


where % = A = Ai = A 2 , and Ac = 

LOpj — ujj. The single-particle part of the effective Hamil¬ 
tonian density has the familiar form of the equal Rashba- 
Dresselhaus SO coupling: 


- 2m L 
+ ^hSax 


Vl - {-idz -b + V;xt(r) 


(4) 


with the effective two-photon detuning and Raman cou¬ 


pling given by (see Appendix B) 


S = - r?) 

A(A2 + ac2)2W2 


Q.R = 


2%2(A2 - «2) 


A(A2-bK2)2 


1 - 


2 % 2 Ac 


A(Ac2-k2) 


mm- 


(5) 


Before proceeding further, consider briefly some realis¬ 
tic order-of-magnitude estimates for various parameters 
used in the theory based on current experiments in ul¬ 
tracold atomic gases and cavity QED. The first exper¬ 
imental realization of a synthetic SO coupling was car¬ 
ried out on ®^Rb atoms using two counter-propagating 
Raman laser beams with wavelength Xr = 804.1 nm 
[Er = 2.33 X 10“^° J) [8]; the two-body interaction 
strengths for the desired pseudospin states of ®’^Rb atoms 
are reported to be gx = 5.009 x 10“®^ Jm^ and 52 = ffi 2 = 
4.986 X 10“®^ Jm3. With typical average BEG densities 
n of order 102° — 102^ m“3 [2], one obtains grfi/ER ~ 1. 

One might reasonably expect interesting physics to 
emerge when the strength of cavity-mediated interac¬ 
tions becomes comparable to the intrinsic inter-particle 
interactions, i.e. when VUr/gr ~ 1- Most experimental 
work is focused on the strong-cavity limit, where ^ ^ k] 
typical atom-cavity coupling and cavity decay rates for 
8'^Rb are ~ %e ~ 10k ~ 27r x 10 MHz [32, 33]. 
One can attain VUr/gr ~ 1 by choosing A ^ 26 THz, 
rji = ri 2 = —Ac = 10 MHz (for example, Ac « 28k and 
rj R3 2.2k in Ref. 32), and a volume V = mm^; for 
these parameters one also obtains MIr/Er ~ 4 x 10“^. 
The weak coupling regime relevant to the present work 
can be attained by increasing the value of k, for ex¬ 
ample by decreasing the reflectivity of the cavity mir¬ 
rors. Choosing k ~ 27r x 100 MHz one can never¬ 
theless ensure VUr/gi ~ 1 choosing a larger volume 
V = 10“^ mm3 g^g stronger pump fields and cav¬ 

ity detuning r]i = r ]2 = —3Ac = 3 GHz; these choices 
yield MIr/Er ~ 4 x 10“2. Further increasing the driv¬ 
ing field intensities up to ?7i = 772 = 15 GHz at the fixed 
Ac = —1 GHz results in cavity-mediated interactions 
that are an order of magnitude larger than the two-body 
interactions VUr/gi ~ 30 while UXIr/Er ~ 1. 

In Appendix B, which discusses the adiabatic elimina¬ 
tion of the cavity fields and the origin of the long-ranged 
cavity interactions, quantities such as %^/[A(Ac -b in)] 
and Nr/[A{Ac -b in)] are assumed to be small. Using 
the weak-coupling values considered above and assuming 
a typical average BEG particle number Nr ~ 10^, it is 
straightforward to verify that both /jA(Ac -biK)] <C 1 
and '^Q Nr/\A{Ac -b fK)] ~ 10“^ ^ 1. Making use of 
K ^ Ac and defining ^ = 2%^/AAc <C 1, one can write 


n 


R ‘ 


Ac ' 


A, 


- vl ); u, 


ds 


n -n ~ 2 

01(2) — 0=p(±) ~ ^ 771(2)- 


he 

Ac 


mm] 


( 6 ) 


If 771 = 772 then hS = 0 and C/ds = C^i( 2 ) = ^=f(±) 
Ui/h^R = ^ ^ 1. Alternatively, if both pump fields are 
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non-zero {rii,r ]2 7 ^ 0), then defining dU = U 2 — Ui one 
obtains 6U/h6 = Uds/Kln = ^ «C 1. These relations will 
be important below when choosing parameters for the 
theoretical calculations. 


III. GROUND STATE AND EXCITATIONS: 

ANALYTICS 

The above analysis indicates that as long as 771 and 
772 are not too different from one another then 5 

in the following we therefore restrict calculations 
to (5 ~ 0. The effective single-particle Hamiltonian 
can be diagonalized, and expressed in the form = 
A=± with single-particle energy dis¬ 

persion relation 


e'±(k) ^ ^ = P + 1 ± ^ 4^2 + (7a) 

and spinor eigenstates 



where ‘-I-’ and ’ designate the upper and lower band, 
respectively, and sin 20k = ^r/\J The unitless 

parameters k = and O/j = hflu/Eji are defined for 

convenience, where En = h^kj^/2m is the recoil energy. 
Recall that using experimentally motivated parameters 
as discussed toward the end of Sec. II, one can choose 
Clfi ^ 0{1). Note that in deriving this result we have 
assumed that the condensate is confined in a box poten¬ 
tial with negligible occupation of transverse momentum 
states, i.e. k = (0,0, kz). In fact, the nature of the trans¬ 
verse confinement is not important in the current work; 
for example, instead assuming a strong radial oscillator 
potential V{p) = muJ^pP^ (2 one would simply replace 
by k1 + hwp/Eji under the assumption that the conden¬ 
sate occupied the ground state of the radial oscillator. 

The energy dispersion with respect to kz consists of 
two bands with a band gap of at the origin k = 0 . 
The lower energy band e_ (k) is a symmetric double well 
along the kz direction with the two minima located at 


kz = ±fco = ± \/l - ^|:/4, (8) 

for kin < 2 , and it has a single minimum at fcz = 0 
when kin > 2 (the minima along the other two direc¬ 
tions always occur at ky = 0). The operators $(k) = 
(<^+(k), (,f>-(k))^ annihilate a boson at momentum k in 
the upper and lower bands and are related to the field op¬ 
erators through 4'(r) = A=± "'‘^A(k)‘^A(k). Note 


that the laboratory-frame bosonic field operators ^'(r) 
(which gives the observable atomic density distribution) 
are related to ^(r) by the unity transformation ^ 2 , be. 

4r(r) = '^2^#(r). 

The single-particle ground state of the symmetric dou¬ 
ble well (i.e. when kin < 2 ) is two-fold degenerate; the 
atom is either in the left minimum at k = —kg = 
( 0 , 0 , —fco) or the right minimum at k = kg = ( 0 , 0 ,/cg). 
The non-interacting IV-particle ground state, when the 
cavity-mediated interactions are also absent, is therefore 
{N + I)-fold degenerate (any number of pseudospin-up 
atoms, up to iV, can reside in the left well). Nonetheless, 
the two-body and cavity-mediated interactions compete 
with each other to lift this degeneracy. 


A. Variational Approach 


In order to determine the nature of the ground state, 
we assume the following ansatz for the BEC condensate 
wavefunction, 


1pl 

_7/>2_ 

= A/f7|cie-*'=«^ 

cos 0 ko 
- sin 0 ko 

-f 020 *'=“^ 

sin 0 ko 
- cos 0 ko 







where fcg = fc/jfco and h = N/V is the average particle 
density, with N and V being the total particle number 
and volume, respectively. The variational parameters are 
Cl and C 2 with the normalization constraint |cip-|-|c 2 p = 
1. Once they are determined, one can find the relevant 
ground-state quantities such as the total density 7i(r) = 
|■i/:l(r)p -I- |7/;2(r)p, and the magnetization per particle 

s^(r) = [IV-iMP - IV’2(r)P]/h: 

n{v) = h [I -I- 2 |ciC 2 | cos(2fcgz -|- 7 ) sin20ko], (10) 


= (|cip - |C2|^) COS20ko, (11) 


where 7 is the relative phase between ci and C 2 . Note 
that the magnetization Sz is homogeneous while the to¬ 
tal density 7i(r) exhibits fringes in the 2 : direction pro¬ 
vided that C 1 C 2 7 ^ 0. Constraining kin < 2, one can write 
sin 20 ko = ^ 2 i ?/2 and cos 20 ko = ^ 0 ; then these take the 


simpler form n{z) = n 


I -I- Oj 7 |ciC 2 | cos( 2 A:o 2 ; -|- 7 ) 


and 


S. = fco( 2 |ci |2 — l). The energy functional £'[ 01 , 02 ] = 
£0 + ^'int is obtained from Eq. (2) by replacing the field 
operators ipr with the corresponding condensate wave- 
functions ^pr■ This yields £g = —iV£flfifj/4 and 
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-E-int, — 


|sgn(5fi)+52 + 4C/i + 2(5t/- 2C/ds^K + 2512 - sgn(5i) - 52 + 4 (^C/^ss - - 2(5i/ - 

1 / ~ \ 1/2 

+ 2 (|cil^ - |c2p) (4 - fj|ej [sgn(5i) - 52 - 2 ( 5 C/] 

- 2|ciC2p sgn( 5 i) + 52 + 4C/i + 2SU - 2gi2 - (3sgn{gi) + 852 + SUi + 4SU - 2gi2^ 


? , ( 12 ) 


where the two-body interaction strengths are rescaled by 
I51I (for example 52 = 92 /\gi\) and the cavity-mediated 
interaction strengths are rescaled by |5i|/V^ (for exam¬ 
ple Ui = VUi/\gi\). In the above equations we have 
defined 2Uss = U± + Uzp and SU = U 2 — Ui, and 
sgn(5i) = 51/I51I = ±1 denotes the sign of 51. Again, 
recall that using experimentally motivated parameters as 
discussed toward the end of Sec. II, t/i(2) ~ ~ 0{1). 

Eq is the single-particle contribution to the energy and 
is independent of c,, as expected. Minimizing Ei^t with 
respect to Ci determines the ground state of the system. 
The parameters Ui and SU (or U 2 ) are the only cavity- 
mediated interaction parameters having an effect on the 
ground state. 

Consider first the simplest case where 52 = sgn(5i) 
and SU = 0, so that only that last line of Eq. (12) con¬ 
tributes to the interaction energy. Then the energy is 
minimized either with (ci,C2) = (1,0) or (0,1), or with 
Cl = C2 = l/-\/2 (neglecting relative phases). The first 
solution set corresponds to all atoms condensing in a sin¬ 
gle minimum of the single-particle energy dispersion (i.e. 
a single plane wave with wave vector —kg or ko), la¬ 
beled the plane wave phase (PWP). In the PWP the to¬ 
tal density is uniform. The magnetization takes the value 
Sz = ±ko = ±(1 — 17^/4)^/^, with the upper (lower) sign 
corresponding to ci = 1 (ci = 0). For small 17^ the 
magnetization approaches unity. Note that the PWP is 
twofold degenerate; that is, all atoms can condense in 
the left (ci = 1) or right minimum (c2 = 1). The sec¬ 
ond solution set corresponds to atoms condensing into 
a superposition state of plane waves. It is characterized 
by the broken translational symmetry and the resulting 
density n(z) = n[l + |l7i^ cos(2fcoz + 7)] exhibits spatial 
variations in the z (i.e. SO-coupling) direction, so this is 
referred to as the stripe phase (SP). In this phase the den¬ 
sity oscillations have greatest contrast for large 17—)■ 2. 
The SP magnetization Sz is zero. 

The SP solution yields a lower energy than the PWP 
solution when term in square brackets in the last line of 
Eq. (12) is positive. (Recall 52 = sgn(5i) and (517 = 0 so 
that the middle line vanishes identically.) The cavity in¬ 
teraction strength that favors the SP solution is therefore 


Ui > Ui^, where 

Tjo ^ 8 [ffi2 - sgn(gi)] - [ffi2 - 3sgn(gi)] 

4(4-i7|) ’ ^ ’ 

is the critical cavity interaction for the SP-PWP tran¬ 
sition. In the limit of small (Ir, this becomes [/)(, ~ 
4[512 - sgn(5i)] + ^[gi 2 + sgn(5i)]f7|.^ If 512 = sgn(5i) 
the SP is favored for any non-zero, positive cavity inter¬ 
action in the limit > 0. In the other hand when 

flu —> 2 and 512 ^ — sgn(5i), the critical cavity interac¬ 
tion Ui^ diverges and SP is only favored for very large 
positive cavity interaction. 

It is important to verify that the total interaction en¬ 
ergy, Eq. (12), remains positive; the system is stable 
only if d'^Eint/dN'^ > 0. Let us examine this first in 
the SP where ci = C2 = l/v^, for a special case where 
Uds = Uss = Ui (and 52 = sgn(5i) and (517 = 0 as before). 
One obtains 

Eint = [512 + sgn(5i)] ^8 + i 7 |j^ 

+ ii7i(2-f7H)'|. (14) 

Surprisingly, the SP is energetically stable for two- 
component attractive BECs in the presence of spin- 
orbit interactions as long as the inter-species interaction 
strength is sufficiently large and positive. Substituting 
the critical cavity interaction [7)(, into Eq. (14) yields the 
constraint 

57/? (2 — 17/?)^ — 12 

In the limit of 17/? —>■ 0, for the lowest possible val¬ 
ues of the cavity interaction favoring the SP phase Ui > 
^ic = 7[ffi2 — sgn(5 i)], the SP is energetically stable 
as long as 512 ^ 0, with no constraint on the sign of 
the intra-species interaction strength. Thus, the infinite- 
range cavity-mediated atom-atom interactions stabilize 
attractive two-component BECs against collapse, even in 
the absence of a confining potential. For larger values of 
Ui even the inter-species interactions can be attractive. 

The coefficient of Ui in Eq. (14) is strictly positive. 
Therefore, for a given parameter set {sgn(5i), 512,57/?} 
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one can choose arbitrary large positive values of the cav¬ 
ity interaction strength to strongly favor SP without com¬ 
promising stability (i.e. to satisfy Ui > 1 /°^, while ensur¬ 
ing that -Eint ^0). In other words, the minimal cavity 
interaction Ui which favors a stable SP satisfies 


Ui > max < — 


[512+ sgn(5i)] (8 + Cl%^ 

4(2-nRy 


,U' 


Ic 


(16) 


The stability of PWP can be investigated in a similar 
manner. The plane wave phase is favored when Ui < !/((.. 
The positivity constraint of the interaction energy in the 
PWP 


Eint — 


N^lgi 

2V 


sgn(gi) -h I [gi 2 - sgn( 5 i)] 


uy 


2-n 


> 0 , 


(17) 


iS" 



FIG. 1: Phase diagrams in the (a) and (b) {Ui,gi 2 } 

parameter planes. The stripe and plane-wave phases are de¬ 
noted by back and white, respectively; dark (light) grey in¬ 
dicates the regions where the SP (PWP) is nnstable. (a) 
Phase diagram for sgn(pi) = 52 = 1 and different values of 
gi 2 = 0.1, 1, and 2. (b) Phase diagram for sgn(gii) = §2 = —1 
and CIr = 0.1. 


imposes a lower bound in the cavity interaction 

- 8^8i.to) + l!;.2-sgn(g.)|fi|, ^ ^ 

8 (2 - 

beyond which PWP is unstable. Thus, even the PWP 
becomes energetically stable for attractive spin-orbit cou¬ 
pled two-component BECs if the cavity-mediated inter¬ 
actions are judiciously chosen. 

Figure 1 depicts the phase diagrams in the {Ui,flR} 
and {Ui,gi 2 } parameter planes. The phase diagrams are 
comprised of two physical regions: the SP and PWP, de¬ 
noted by black and white in Fig. 1, respectively. The 
dark (light) grey indicates the regions where the SP 
(PWP) is energetically unstable. Figure 1(a) shows 
the phase diagram in the parameter space for 

sgn(gi) = ^2 = 1 and different values of gi 2 . The stripe 
phase is favored over an ever-larger parameter space as 
Ui increases as long as |Ofl| < 2 to assure the existence 
of a double-well single-particle dispersion. This general 
trend is also evident from Fig. 1(b), the phase diagram 
in the {Ui,gi 2 } parameter plane for sgn((/i) = 32 = —1 
and constant flu = 0.1, where Eq. (13) reveals that the 
phase boundary is linear in gi 2 for fixed Qr. 

Relaxing the constraint considered above that 6U = 0 
in Eq. (12), one can prepare any arbitrary superposition 
state, i.e. arbitrary ci and C 2 subject to |cip -I- |c 2 p = 1. 
The plane-wave phase is no longer degenerate; rather, the 
minimum favored depends on the sign of SU. Figure 2 
shows the dependence of |cip in the {C/ 12 , ^r} plane for 
sgn(gi) = g 2 = 6U = 1, and gi 2 = 2. Under these condi- 
tons the SP with |ci| = |c 2 | is found only for very large 
Ui ^ Uic, i.e. far from the SP-PWP phase boundary 
Eic. Whereas for Ei —>■ C/{{, |ci| increases monotonically 
until the PWP with |cip = 1 is attained for C7i < Uic 
(note that the critical value Eic — Ui^ and is weakly 
dependent on 6IJ, as discussed below). The plane-wave 


phase begins to be unstable in the left bottom corner of 
this figure. 

The magnetization Sz = fco (2|cip — l) as a function 
of C/i is illustrated with the black solid curve in Fig. 3 
for sgn( 5 i) = g 2 = SU = 1, gi 2 = 2, and (Ir = 0.1. For 
contrast, the magnetization when SU = 0 is also shown 
(blue dashed curve). Note that while the sign of the 
magnetization in the PWP is arbitrary for the SU = 0 
case (a spontaneously broken symmetry in the ground 
state), in the present case the sign of Sz always follows 
that of SU. On the PWP side, the magnetization is fixed 
at its maximal value Sz = fco! for C/i > Uic on the SP 
side, the magnetization decreases sharply before reaching 
an asymptotic value deep within the SP phase. 

For small SU and ^Ir, the SP-PWP phase transition 
occurs at almost the same value of the critical cavity 
interaction C/^^, = 0.5 obtained using Eq. (13) which as¬ 
sumed SU = 0. Near the phase transition point on the SP 
side, one can write ci = 1 — and C 2 = y/2x, where x <C 
1 and cf -I- C 2 — 1 -I- 0{x‘^). Setting the term proportional 
to in Ei„t[ci = 1 , C 2 = 0 ] - Ei„t[ci = 1 - x^, C 2 = v^x] 
equal to zero yields a modified critical cavity interaction 


Uic = C/i°e 


2-(4-U|j)1/2 




2 

R 


SU. 


(19) 


In the small Hr limit this may be simplified to Uic — 
5 [512 - sgn( 5 fi)] -k rU + sgn( 5 i) -k SU]H% which is the 
same critical cavity interaction [/{{ obtained above in the 
small Hr limit, save for the 5C/-dependent correction. 

The behavior of the magnetization for Ui > Uic sug¬ 
gests that one can define the order parameter for the 
stripe phase to be E = 1 — Sz/ko = 2(1 — c{). As desired, 
this vanishes in the PWP (here we only consider a PWP 
with momentum —kg) and takes a nonzero value in SP. 
The order parameter is shown in the inset of Fig. 3. The 





7 



0 1 2 


Ci/f 

FIG. 2: (Color online) Density plot of |cip in the 
parameter plane for sgn(gi) = g 2 = SU = 1, and gi 2 = 2. 
The plane-wave phase begins to be unstable in the left bottom 
corner. 


discontinuity in the derivative of P with Ui suggests that 
the SP-PWP quantum (zero-temperature) phase transi¬ 
tion is second order. It is therefore of interest to deter¬ 
mine the (mean-field) exponent f3 for the order parame¬ 
ter P in the vicinity of the transition point. Substituting 
Ui = Uic + X info the energy functional Emt and mini¬ 
mizing it with respect to Ci yields 


2SU 

(4-^1 

vl/2 

+ x{ 

^4 ^1;) 

2iU{ 

^4 — 

1/2 

+ 2x 

(4 - 


( 20 ) 


The order parameter P = 2(1 — cf) computed using this 
expression for Ci is illustrated as the green dashed curve 
in the the inset of Fig. 3, and is in excellent agreement 
with the numerical results of the variational approach, 
shown as the black solid curve. Taylor expanding ci in 
Eq. (20) for small y and 17up to first and second order, 
respectively, one obtains ~ 1 — x/2i5[/ (the term pro¬ 
portional to X^ji is also omitted). This yields the mean- 
field order parameter Pmf = 2x/(5C/ = 2(Pi — Uic)^/SU 
and a critical exponent /3 = 1. The behavior of the order 
parameter near the transition point fits well to P, as is 
shown by the orange dashed curve in the inset of Fig. 3. 

In principle, it is not valid to consider SU ^ 0 while at 
the same time assuming that S = hS/En = 0. Rather, 
a rji ^ r ]2 0 but rji ^ 772 , then Eqs. (6) state that 

S ~ SU whenever Ui ~ That said, in Fig. 3 the 
parameters are chosen so that CIr = 0.1 ^ SIJ = 1. One 
can therefore expect <5 ^ SU by a similar ratio, which 
again justifies neglecting it. 

Consider briefly the effect of keeping a non-zero but 
small value of S. The single-particle dispersions of the 





FIG. 3: (Golor online) The magnetization as a function of 
Ui shown as the black solid curve for sgn( 5 ri) = §2 = SU — 1, 
gi 2 = 2, and CIr = 0.1. The blue dashed curve represents the 
magnetization when SU = 0. The red dotted curves are the 
magnetization computed from solutions of the coupled Gross- 
Pitaevskii equations in the SP and PWP assuming Uss = 
Uda = Ui for the same parameters as the solid black curve, 
and \gi\n/Eii — 1. Inset: the SP order parameter P is shown 
as a function oiUi (black curve); an analytical approximation 
(dashed green curve) and the behavior near the critical point 
(orange dashed curve) are shown for comparison. 


spin-orbit Hamiltonian (4) become 


e±(k) 



( 21 ) 


rather than the expressions given in Eq. (7a). The asso¬ 
ciated (orthogonal) eigenvectors have the same form as 


Eqs. (7b) but now sin2dk = ^r/\Iz [^^z + S] + ^r- 


'V^ 


For d 7 ^ 0, the lower double-well dispersion curve e_ is no 
longer symmetric; rather, the right well is lower (higher) 
when (5 > 0 (d < 0). Thus, in the absence of particle in¬ 
teractions a PWP is energetically favored in one well or 
the other with no ambiguity. The presence of S precludes 
a simple form like Eq. (8) for the location of the energy 
minima, but in the limit when both ^Ir <C 1 and d <C 1 
one obtains 



The lowest-order contribution of d is a correction to the 
coefficient of the already small H^^-dependent term, and 
therefore the value of fcg is well-approximated by assum¬ 
ing 5 = 0. Likewise, the BEG approximation consists 
of kz with fco; because ikz —>■ 4fco « 4 3> 5 in the ex¬ 
pressions for the single-particle energies and eigenvectors 
above, S can be similarly neglected in the calculations. 













B. Coupled Gross-Pitaevskii Equations 


real, Ui = Uss = C4 js, and sgn( 5 i) = g 2 , one obtains 


While the variational calculation discussed in the pre¬ 
vious section has revealed that a ring cavity can stabi¬ 
lize stripe phases in interacting spin-orbit coupled Bose- 
Einstein condensates, it is important to verify the re¬ 
sults using a more rigorous approach. In this section, the 
coupled Gross-Pitaevskii (GP) equations are derived for 
both PWP and SP ansatze and the ground state proper¬ 
ties are obtained from their solutions. 


1. Plane wave phase 


\9i\ 

Er 


{gi 2 - sgn(gi)) (V-a 


^l) ^ 1^2 + t/i 


± 4fcoV'lV’2 + (V '2 - V'l) = 0. 


( 25 ) 


For the plane-wave phase, both 'ipi and ip 2 are assumed 
to be constant, so that ijji +' 4)2 = n and 'tpl — = fiSz- 

Inserting these into Eq. (25) gives 




+^ke) 


Sz^^\gi2 - sgn( 5 ] 


Er 


+ 2sz I + f7] 


Iff 


i\n 


Er 


= 0 . 


(26) 


The GP equations can be obtained directly from the 
many-particle Hamiltonian (2): 


n —+ S'llV’iP + ffi2|'*/’2p + 2UiNi + 2UdsSx V’l 
2m 

+ [H^Ir, + 2UssS- + UdsEf'\'4’2 = mV’i, 

-—Aa -b <72|V’2p + ffi2|V’iP + 2C/aAa -b 2UdsSx V’2 
2m 

+ \hflR + 2UssSd- + UdsN'l'P’i = mV’27 (23) 


where Ai = —-b {—idz + kR^ and Aa = —-b 
{—idz — kR)^ and the BEG wavefunctions for the two spin 
components are denoted by V’i( 2 ) rather than V’i( 2 )(r) 
to save space. These equations can be simplified in the 
plane-wave phase by assuming homogeneous wavefunc¬ 
tions V’T(r) = where the upper (lower) sign 

corresponds to a condensate in the right (left) minimum. 
The GP equations are then recast as 


/i- (fco± 1)^ 
^ \9ii 

■ Er 


V'l - ^Rtp2 

(sgn{gi) + 2Ui^ + (ffi 2 + 2Ass) |V ’2 

-b Ads (2|AiP +1^2^) V’2 + AdsAiA^ 

V'2 - ^Rlpl 


V’l 


fl- (fcoT 1)^ 

\l 

■ Er 


(^ff2 + 2Aa^ IV'2|^ + ^5i2 + 2Ass^ IV’il 


V’2 


+ Kis (IV’il +2|V’2| ) V’l + KisV’2V’l 


(24) 


where again the upper (lower) sign in each equation cor¬ 
responds to a condensate in the right (left) minimum, 
and the chemical potential is expressed in recoil energy 
units, p, = p/Er. 

The chemical potential can be obtained from the first 
of Eqs. (24) and then substituted into the second. Under 
the assumption that both condensate wavefunctions are 


When Ui = 0 and ^r ss 0, this expression is approx¬ 
imately correct when « 1, consistent with the vari¬ 
ational results in this regime. Recall that in the varia¬ 
tional approach, the magnetization Sz = fco is constant 
[c.f. Eq. (11)], solely determined by CIr. Unlike the vari¬ 
ational result, however, it is immediately apparent from 
the second term in Eq. (26) that the magnetization must 
decrease monotonically as Ui is increased. 

The magnetization Sz obtained via numerical solution 
of Eq. (26) is shown as the red dotted curve in Fig. 3 for a 
condensate in the left well (i.e. choosing the lower sign) of 
the PWP for Ai < Aic. Parameters are Ui = Ugs = Ads, 
sgn(ffi) = 92 = \gi\n/ER = SU = I, gi2 = 2, and 
(Ir = 0.1. As expected, the magnetization decreases 
monotonically with Ui from its maximum at Ui = 0. 
The difference between the results of the two methods 
has its origins in the fact that the variational ansatz, 
Eq. (9), is a single-particle wavefunction which satisfies 
the GP equations in PWP only when all the two-body 
and cavity-mediated interactions are zero. In principle, 
the variational ansatz could be remedied by allowing both 
ko and dko to be variational parameters [12]. The depen¬ 
dence of the solution of GP equations on the two-body 
and cavity-mediated interactions will be investigated fur¬ 
ther in Sec. Ill C 1 in the calculation of the elementary 
excitations in the PWP. 


2. Stripe phase 

The momentum dependence of the condensate in the 
SP is not as readily apparent as it is for the PWP. It 
is therefore convenient to instead construct an effective 
low energy Hamiltonian by first mapping the complete 
Hamiltonian (2) into the lower band and then deriving 
the low energy coupled GP equations [14, 34]. This is 
reasonable because the occupation of the upper band 
e+(k) can be assumed to be small at low temperatures 
/cbT <C EPIr. Furthermore, only states in the vicinity of 
the two minima ikg will be occupied. 

The field operators ’®'(r) can then be expanded in the 
lower band basis around the two minima (recall that 
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4 >_{k.) is the two-component spinor in the lower band): 


q<qc 


+ (ko + q)^_ (ko + q)l, (27) 


dispersion ee(k)/iili^ = k‘^-\- — (2^/4)fc^. It is impor¬ 

tant to note that the lowest single-particle energy state 
for both dressed psendospins is the k = 0 momentum 
state, not k = ±ko as it was for the actual pseudospins. 
Then the effective low energy GP equations for the SP 
can be obtained from iJg, Eq. (30): 


where the sum over q need only be taken up to some 
maximum qg. Approximating the spinor ^_(±ko -fq) ~ 
(p_ (±ko) in the limit flu 2 and defining the new oper¬ 
ators (,3i/(q) = i^_(-ko-bq) and (^ 2 '(q) = i^-(ko-fq) [14], 
the field operators read 

4'(r) = e"*‘"“''’</>_(-ko)V:i/(r) -be*’^‘’''’^_(ko)Vi 2 '(r), 

(28) 


where (q)- In the small flji limit 

and keeping terms only up to second order in flu and 
noting that ko — {1 — f 2 ^/ 8 )fc_R, the field operators can 
be further simplihed to 


'tPiir) 


.^ 2 (r)_ 



\„-iknz 
32 m 

Or -ikoz 
4 ^ 


-(1 


„iknz 
4 e 

_ Ha'ipifcoz 
32 m 


^i'(r) 

> 2 '(r). ’ 
(29) 


Note that the lab-frame pseudospin field operator ipT 
maps correctly to the corresponding dressed pseudospin 
field operator ipT' in the CIr —> 0 limit; recall that 

’®'(r) = '^ 2 ^’®'(r). Substituting Eq. (29) back into the 
original Hamiltonian (2) and only keeping terms to sec¬ 
ond order in flR yields the effective low-energy Hamilto¬ 
nian: 

= / d^r -b -b + g'i2nvn2^ 

+ \u[Nl + + U[^NyN2., (30) 

where 4^ (r) = (■i/)i/(r),' 02 '(r))''', as before Nt' = 
f nTf(r)d^r = f ipl,(r)ipTi(r)d^r is the total atomic num¬ 
ber operator for the dressed pseudospin r' € {!', 2'}, and 
we have introduced the dressed interaction parameters 

9t = 9t't' = 9t ~ ~^^9t ~ 9l2)^R^ 

9i2 = 9i' 2' = 9i2 + g(5i + 92 )^% 

C/; = Ur'r' = 2Ur - Uos^R - \{Ur - Uss)^^, 

1/(2 = Uv2> = -Uds^R + ^(C/i + 1/2 + 2Uss)fl% (31) 

with r e {1, 2} and t' € {!', 2'}. 

The single-particle part of the effective low energy 
Hamiltonian = (—li^/2m)[V^-b(l —H|j/4)i9^] can be 
easily diagonalized [14], yielding the effective low energy 


\ 9 i + Ui) 

li’i'f + 

(912 + H12) 

1 IV' 2 '[^ 

Ipl' 

= mV'i' 

(92 + uCj 

[V' 2 '[^ + 

{912 + H12) 

IIV'i'P 

4 ’ 2 ' 

= MV' 2 ' 


(32) 


where the dressed pseudospin wavefunctions ipT' are as¬ 
sumed to be homogeneous and unitless parameters have 
been introduced for convenience: g'^ = 9 '^/\gi\, §[2 = 
9j2/\9i\, K = Et/)./[gi[, and 1/(2 = VU[ 2 /\ 9 i\- Here 
jl = fJ-/\gi\ which has units of inverse volume. These 
algebraic equations have the solution 

2t/2 + 52-5i2-i (H 1 + 3 C/ 2 ) 

nv = - jz - ' N / \ 

9i + 92~ “2912 + 2 \JJi-\- 1/2 j ^1 — jfljij 

2Ui +91— 9 i2 — I ( 3 H 1 + C/ 2 ) 

” 2 ' = —:-:- r- —^tt— 

9i + 92~ 2ffi2 + 2 (jJi + C/ 2 ) ^1 — 4 HI.) 

where ny + n 2 ' = n. Note that although the GP equa¬ 
tions for the SP, Eq. (32), depend on the cavity parame¬ 
ters Uss and C/ds, these solutions do not; rather, Ui and C /2 
are the only cavity interaction parameters that affect ipr’, 
consistent with the variational approach of Sec. HI A. 

The dressed magnetization si, = (ni' — n 2 ' )ln can eas¬ 
ily be obtained from Eq. (33), and the actual magneti¬ 
zation Sz = s(,(l — H|j/8) up to 0(H|j) is found using 
Eq. (29): 


Sz = 


~g'2-9[+26ij(l-lfll)] (l-|H|) 


g'l +92~ 2g'i2 + 2 (Ui + t/ 2 ) ^1 — qHfl) 


(34) 


The SP magnetization Sz is displayed as a function 
of t/i(^ Uic) in Fig. 3 with the red dotted curve for 
sgn(gi) = 92 = SU = 1, 912 = 2, and ^r = 0.1. The 
behavior is indistinguishable from the magnetization ob¬ 
tained from the variational approach, Eq. (11). The crit¬ 
ical cavity interaction for the SP-PWP phase transition 
can be obtained from Eq. (33) by setting nR = n (or 
setting si, = 1 ): 


f+ — 

'^ic — 


4(4-H2j) 


912 - sgn(gi) - 292 - SU 


+ 8 [gi 2 -sgn(gi)] 


fl 


R 


(35) 


for a phase transition from SP to a PWP at the left min¬ 
imum. Instead setting nR = 0 (or s^ = —1) for a phase 



10 


transition from SP to a PWP at the right minimum, one 
obtains 


77 R _ 


4(4 


gi2 - sgn(5i) - 2g2 - 3SU\ 


R 


gi2-sgii{gi)-2SU\ }. (36) 


Note that when sgn((/i) = g 2 and 6U = 0, the two critical 
cavity interactions and C/^ become equal to the value 
[/((. found using the variational approach, Eq. (13). 


C. Elementary Excitations: Bogoliubov theory 

Thus far we have treated the bosons as classical fields, 
having replaced the field operators with their expecta¬ 
tion values t ' 0 T = In this section, we consider 

the quantum fluctuations of the fields and obtain the 
elementary excitation spectrum using Bogoliubov the¬ 
ory. This is accomplished by writing the field operators 
as ijjr = ’’pT + 5'4 ’t, where Sipr is the quantum fluctu¬ 
ation operator. These expressions are substituted into 
the time-dependent GP equations and the resulting equa¬ 
tions are linearized, i.e. terms are retained only up to first 
order in the fluctuations. One then obtains a set of time- 
dependent coupled equations for S'tpr which yields, after 
diagonalization, the elementary excitation spectrum. 


1. Plane wave phase 


Following the approach taken in Sec. Ill B 1 for the 
PWP, it is reasonable to define the bosonic field operator 


V't (r,t) = Ur + dV’r (r, i) 


(37) 


where ipr are the time-independent, homogeneous solu¬ 
tions of the coupled GP equations (24) in the PWP. To 
consider time-dependent fluctuations around the equi¬ 
librium solutions it is convenient to replace the chem¬ 
ical potential (which is the eigenvalue of the time- 
independent GP equations) by a time-dependent oper¬ 
ator, /i —> /r -f ihdt . The time-dependent fluctuations can 
then be expressed using the usual Bogoliubov approach 
in terms of particle and hole excitations with amplitudes 
and u* respectively. 

Consider t_he specific case of a condensate in the left 
minimum —kg of the double-well single-particle disper¬ 
sion relation; for condensation in the right well one need 
only replace ko in what follows with — fco- Substitut¬ 
ing Eq. (37) into the time-dependent GP equations and 
keeping only linear terms in the fluctuations, one obtains 
the following non-Hermitian eigenvalue equation for each 
value of q: 


/ Mn 

-gii’f 

gi2'P’*i'>p2 + kills 
\ -gi2'iPt'4’2 


gi'ipl gi2'4’i'4>i + kites 

-M 22 -gi2'4>lPpl 

512 ^ 1^2 M33 

-512^1'*/'$ - kites -g2'4’l‘^ 


where 


Mii /22 — Er 9^ T 2(fco ~ l)9z + gi\ipi\'^ ~ kites^^, 

L J iPi 

11 ^ 33/44 = -S'i? T 2 (fco + l)9z + g 2 \i^ 2 \^ — kills'^, 


hiteS = kiln + \gi\Udsn + 2\gi\Ueeipi'ipl. 


(39) 


In deriving the Bogoliubov Hamiltonian (38), we made 
use of the fact that Nr = J ipl(r,t)i’T(i^,t)dr = 
J Itprl'^dr = V\'(pr\‘^ = Nr, becausc Ppr is homogeneous 
by assumption and / 6iprir,t)dr = 0 because the spatial 
integral of either Bogoliubov amplitude or 

h* is zero for any q 0. A similar argument 

ensures that S+ = S+ and S- = S- as well. Note also 
that the chemical potential in Eq. (38) has been elimi¬ 
nated using the coupled GP equations (24). 


gi2'4’i'4’2 

-gi2'ip*i'4’2 - kills 
52^1 

— ik /44 


\ 



fupq\ 


^l.q 

= k(jj{q) 

hl,q 


M2.q 


■U2,q 

/ 

\^2,qy^ 


\^2.q/ 


(38) 


Diagonalizing Eq. (38) yields the spectrum a;^'^(q) of 
collective excitations. The results are shown in Fig. 4(a) 
for the parameters sgn((/i) = 52 = \gi\k/Eii = 1, 
gi 2 = 2, and itn = O.I, when all the cavity-mediated 
interaction terms are zero (Ui = U 2 = Uss = Uss = 0), 
i.e. the system is deep in the PWP. The lower curve ex¬ 
hibits the usual superfluid sound-like linear dispersion 
around the origin = qz/kfi = 0 (around the left min¬ 
imum of the single-particle energy dispersion where all 
the atoms are condensed) and a roton-type minimum 
around (jz — 2. As the cavity interactions are increased, 
the energy of the roton minimum lowers. For parame¬ 
ters Ui = 0.5, 61J = 1.5, Uss = Uds = 0, and the other 
parameters same as in Fig. 4(a), this minimum coincides 
with zero energy (i.e. the excitation energy at the ori¬ 
gin (jz = 0); see the black solid curve in Fig. 4(b). The 
red dashed-dotted curve represents the elementary exci¬ 
tation spectrum for the same values of C/i and 6U but 
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riaF^(<i)IER 



(a) 





(b) 


FIG. 4: (Color online) Elementary excitation spectrum in the 
PWP for sgn(yi) = P 2 = |sri|h/j5_R = 1, gu = 2, and = 
0.1. {U\,U 2 ,Uaa,Uia) = (0,0,0, 0) in (a), and (0.5,1.5, 0,0) 
and (0.5,1.5,0.5, 0.5) in (b) for the black solid and red dashed- 
dotted curves, respectively. 


for ilss = Ud& = 0.5. In this case, Hiles/En [cf. Eq. (39)] 
is somewhat bigger than the bare = 0.1 for the black 
solid curve, so the roton minimum lies somewhat above 
that of the black solid curve around ~ 2 . 

The energy of the roton minimum near ~ 2kR can 
be reduced below zero by further increasing the cavity 
interaction strength Ui. This signals a dynamic instabil¬ 
ity toward the formation of the SP; recall from Eq. (10) 
that the density modulation in the SP has wave vector 
2fco ~ 2/c/j for (Ir —>■ 0. The critical cavity interactions 
for the black solid and the red dashed-dotted excitation 
spectra in Pig. 4(b) are t/ic — 0.5 and 0.53, respectively. 


and these are in good agreement with that of the varia¬ 
tional approach, where Eq. (19) predicts a phase transi¬ 
tion between the PWP and the SP at the critical value 
Uic ~ 0.5 for the parameters sgn(gi) = g 2 = 6IJ = 1, 
gi 2 = 2, and Qr = 0.1 (cf. also Fig. 3). 

If one hypothetically sets t/gg = Uds = 0 in the PWP, 
then the critical cavity interaction Uic obtained from the 
analysis of the elementary excitations and the variational 
method would match exactly with each other for any 
range of parameters. Nevertheless, they begin to devi¬ 
ate from one another as C/gg and Uds become larger and 
larger, because Eq. (19) is independent of these cavity 
interaction parameters while both the coupled GP equa¬ 
tions and the Bogoliubov Hamiltonian depend explicitly 
on them (the latter through hflgg)- That said, we have 
compared the critical phase transition point Uic obtained 
from both the variational approach and the elementary 
excitation spectrum in the PWP and have found that 
when Ui = Uss = Uds they agree with one another within 
a ~ 8 % error for 512 in the range of ~ 0 — 8 , assuming 
sgn(gi) =92 = \gi\n/ER = SU = 1 and CIr = 0 . 1 . 


2. Stripe phase 

The derivation of the Bogoliubov excitation spectrum 
begins with the corresponding time-dependent, effective 
low energy GP equations in the SP [c.f. Eq. (32)]: 


+ + U[Ni, + U[2N2- - 9 )Ur, 

zS^^2' +521^2'!' + C/ 2 ^ 2 ' + Ui^Np - g)U2'- (40) 


As in the PWP case, the low energy field operators are replaced with Ut' (r, t) = Ur + SUt' (r, t) in these equations. 
Here Ur are the time-independent, homogeneous solutions of the effective low energy GP equations in the SP, Eq. (33), 
and SUr (r, t) are the quantum fluctuations. Linearizing Eq. (40) yields the Bogoliubov Hamiltonian: 


/ee{q) + g'llUi'l^ 

g'lUl 

9'l2UvUl' 

g'i2UrU2' 

\ 

/wiCqX 




-9'iU*i? 

9'i2UvU2' 

-e-e{q) - g'i\Ur? 
g'i2UvU2' 

-9'i2UvU*2' 
ee(q) +52lV'2'P 

-g'i2UvU2' 

92 UI 


^I'.q 

U2',q 

= Huj{q) 

Vl',q 

U2',q 

(41) 

\ -9'i2UvU2' 

-9'i2UrU*2' 

-92U2^ 

-ee(q) -52lV'2' 

IV 

\'>J2rq) 


\V2',q) 



which can be diagonalized to give the spectrum of the with 
elementary excitations: 

El = g'lUv + 52 ’^ 2 ', 

D2 = {9W2 - 9?2)'nrn2,. (43) 


(q) = J ee(q) + ee(q) {di ± ^ (42) 


We have again used the fact that Nt' = Nr' ■ 

Surprisingly, the Bogoliubov Hamiltonian in the SP 
does not depend explicitly on the cavity parameters and 
the form of the excitation spectrum coincides with the 
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FIG. 5: (Color online) The speed of sound in the transverse 
direction is shown as a function of U 2 for Ui = 1/4 
(solid curves) and Ui =5/2 (dashed curves). For all curves: 
CIr = 0.4, sgn( 5 i) =^2 = 1, 512 = 0.7, gin/Eu = 1, and m is 
the mass of ®^Rb atom. The insets show the results closer to 
the origin. 


quasiparticle spectrum of a Raman-induced stripe phase 
BEG [14]. That said, the excitation spectrum implicitly 
depends on the cavity parameters t// through nr', as can 
be seen in Eq. (33). Both a;|.^ (q) are gapless and exhibit 
linear dispersion at long wavelengths, the characteristic 
of superfluidity in this phase; the slope of the dispersion 
relation at long wavelength corresponds to the speed of 
sound in the medium. In the transverse direction, one 
obtains 


.,(±) 


dcofiq) 

dq± 


q—>0 




-4i^2, 


(44) 


and the speed of sound in the z (SO-coupling) direction 
is nearly the same for small Cl, y^l — f4|./4. 

Eigure 5 depicts as a function of U 2 for 17i = 1/4 
(solid curves) and Ui = 5/2 (dashed curves), with the 
other parameters fixed to sgn( 5 i) = 32 = 1) ffi 2 = 0.7, 
CIr = 0.4, and gin/Eji = 1. The mass is assumed to be 
that of ®’^Rb. As U 2 is increased above zero, the speed 
of sound in the positive branch (the blue curves) 
first decreases quickly and reaches a minimum around 
6U = U 2 — U 1 ^ 0 for both curves, and then gradually ap¬ 
proaches its asymptotic value. The speed of sound in the 
negative branch (black curves) has the opposite be¬ 
havior, first increasing sharply to a maximum again near 
SU ^ 0 for both curves, before asymptotically approach¬ 
ing zero. The insets show the behaviour of close 
to the origin. The asymptotic behaviour of the speed of 
sound can be understood by noting that for large positive 
U 2 S> Ui, ni> approaches n and n 2 > approaches zero [c.f. 
Eqs. (33)]. As a consequence ZI 2 —t 0 and —>■ 0 while 

—)■ ^g'^fijm. For the solid curves (where t7i = 1 /4), 
the speed of sound in the negative branch ^ becomes 


zero at U 2 — 37, consistent with the value at which the 
dressed magnetization s(, becomes unity for this choice 
of parameters. This signifies an instability toward the 
formation of a different phase. 

The condition that the speed of sound must be non¬ 
negative imposes the constraint D 2 ^ 0. This condition 
marks the onset of a phase transition at the critical point 
9vi = \/ff 152 > which does not depend on any cavity- 
mediated interaction parameters and is solely determined 
by the two-body interactions and This critical point 
is not consistent with the previous results obtained from 
the variational approach, the effective low-energy GP 
equations in the SP, or the elementary excitations in 
the PW which all consistently predict a critical point 
for the PWP-SP phase transition that depends on the 
cavity-mediated interaction parameters. To verify that 
there was not an error in the calculations, the elementary 
excitations were computed directly in momentum space 
by Fourier transferring the effective low-energy Hamilto¬ 
nian (30), and treating the fluctuations around the con¬ 
densate — 0) to second order in i,3i-'(q) for small 

momenta q. The results were identical with the real- 
space analysis, Eq. (42). Interestingly, the critical inter¬ 
species interaction g^!^^ above defines a phase boundary 
between the stripe phase and a phase-separated state 
in Raman-induced spin-orbit coupled BECs [14]. It is 
therefore conceivable that there is another phase between 
the SP and the PWP induced by the cavity interactions, 
whose signature is the observed inconsistency in the crit¬ 
ical point. 


IV. DISCUSSION AND CONCLUSIONS 

In this work we have shown that cavity-mediated long- 
ranged interactions between atoms can profoundly al¬ 
ter the nature of the ground state and the elemen¬ 
tary excitations of a cavity-induced spin-orbit-coupled 
two-component BEG, for ring-type cavities in the weak- 
coupling regime. Specifically, experimentally tunable 
cavity-mediated interactions compete with the standard 
two-body interactions to yield both plane-wave and 
stripe phase ground states. Indeed, positive long-range 
cavity interactions can stabilize fully attractive BEGs 
(condensates where intra-species collisional interactions 
are negative, independent of the sign of the inter-species 
interaction) against collapse in the stripe phase. The col¬ 
lective excitations of the plane-wave phase ground states 
are found to have a distinctive roton-type excitation spec¬ 
trum reminiscent of that of superfluid ^He, which can be 
used as a signature of the phase. The stripe phase has a 
standard linear dispersion relation; the associated speed 
of sound is found to go negative at a critical value of 
the cavity interaction strength, signalling an instability 
toward another (likely phase-separated) phase. The re¬ 
sults suggest that cavity QED, even in the weak-coupling 
regime, can yield interesting new physics for spin-orbit 
coupled BEGs. 
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The results raise interesting avenues for future inves¬ 
tigations. This work assumed a fictional experimental 
configuration where the momentum is a good quantum 
number in the direction of the applied spin-orbit inter¬ 
actions. In reality the condensate would be confined in 
this direction, and even a weak harmonic potential could 
change the physics. While the stripe phase would likely 
remain robust, as it is essentially a weak standing wave 
superimposed on the background condensate density pro¬ 
file, the plane-wave phase has no analog in a conhned 
geometry. Another loose end is the nature of the phase 
hinted at in the limit of a large difference SU between 
the cavity-mediated interactions between the two kinds 
of spin components Ui and U 2 - For large 6U, the sound 
velocity in the stripe phase was found to go negative, a 
signature of the dynamic instability of the phase. While 
other work suggests that this signals a 

However, a few intriguing issues and questions remain 
unclear and deserve further investigations. These in¬ 
clude the inconsistency in the critical phase transition 
point, how the combined SO coupling effect, the two- 
body interactions, and the cavity-mediated long-ranged 


interactions change the superfiuid-Mott-insulator phase 
transition as well as the nature of magnetic orders in 
the Mott-insulating regime when an optical lattice im¬ 
posed inside the cavity. Furthermore, whether it is pos¬ 
sible to have a superfluid-Mott-insulator-like phase tran¬ 
sition with solely the cavity-mediated long-range interac¬ 
tions, whether there is more interesting physics in strong¬ 
coupling regime, and how the cavity fields are affected by 
the atoms. Some of these questions are the subject of our 
current works with some promising preliminary results 
and will be published elsewhere. 
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Appendix A: Adiabatic Elimination of the Atomic Excited State 


We first express the single-particle Hamiltonian density Eq. (1), in the rotating frame of pump lasers [29] by 
applying the unitary transformation 


= exp |z (^A\Ai - (Taa) Wpi + {a\A 2 - abb^ a;p 2 j t|, 


to obtain 




rp^ 

2m 


+ Fext(r) 


M' \ ^ / A A \ 

13x3 H-^ yo'aa ~ CTbb) — — (Aal -f Aa2j O'ee + U 


— /i(AeiA|Ai -|- AC 2 A 2 A 2 ) -f ih T — H.c. 


+^&ee-*'=«"i2ae&) +H.C. 

(Al) 


where we have defined the atomic and the two-photon (or relative-atomic) detunings 

Aal — ^pl ^a 2 — ^p 2 ^ — (^pl ^p 2 ) ^ — ^al Aa 2 j (A2a) 

and cavity detunings 


Acj - ^pj ^jl j - ^5 


(A2b) 


with respect to the pump lasers. Let us now assume that the detunings Ai = wi — £ea/fi = —Ad -f Aai and 
A 2 = W 2 — Seb/fi- = —Ac 2 + Aa 2 are large compared to Sba/fi = (sb — £ 0 )/^ so that we can adiabatically eliminate 
the dynamic of the atomic excited state |e) from the Hamiltonian (Al) and obtain an effective Hamiltonian for the 
ground pseudospins {1, 2} = {&, a}. Hollowing the standard adiabatic elimination procedure [23, 35], we first hnd the 
Heisenberg equations of motion ihoer = [oeT,'H!^^'^\ for dea and deb, and then (after transferring to slowly rotating 
variables) set them equal to zero to find the steady-state solutions oi‘^^ and . After substituting these steady-state 
solutions back in (Al) and dropping terms diagonal in we arrive at the single-particle Hamiltonian density 

for pseudospins 


^SO — 


r „2 




I + eiCTii +£ 2*^22 A\A2cr2i^ 


tt' 

-^cav’ 


(A3) 
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where 


-^cav ~ + ^02^2^2) + 


(viA 


ri 2 AX] - H.c. 


and 


£l = - 


hS' 

~Y 


A 2 


{A\A2 + 2 )’ 



/ it i 


(A4) 


Here, '^ae^be is the two-photon Rabi frequency and I = 12 x 2 is the identity matrix in the pseudospin 

space. Note the hat on Sr, implying that it depends on the cavity field operators. After transferring to the co-moving 
frame of the cavity modes by applying the unitary transformation ^2 = e~'‘^^Aii-<^ 22 )z |.]^g Hamiltonian density 

(A3), we obtain the SO-coupled single-particle Hamiltonian density 


^SO 


1 

2m 


+ [pj, - hkR{<J22 - CTll)]^! + Rext(r)II - 


E 

r=l,2 


e-rGTT + 


fl( 


AIAiCTi2 


+ 41 


20'21 


H' 

■^cav 


(A5) 


One can identify hkR{a 22 — o'ii) with eAl/c as in the minimal coupling Hamiltonian, that is, eA^jc = hkii{a 22 — o’ii) = 
—hk^az, where Uz = un — CT 22 is the third Pauli matrix. Nonetheless, we emphasis that here A* is a matrix acting in 
the internal pseudospin states, in contrast to the ordinary vector potential whose components are scaler fields. Then 
the single-particle Hamiltonian reads 

Hso^ = ^ + Kxt(r)l] irNz + hn'j,(^AlAiS+ + 44^-) + H'av, (A6) 

where ^(r) = (^/;i(r), 4 (r))^ are the bosonic field operators, 4 = /is the total atomic number 
operator for pseudospin t, N = Ni + N 2 is the total atomic number operator, and 4 = 4= / (r)4(r)d^r’ are 

the collective pseudospin raising and lowering operators. 


Appendix B: Adiabatic Elimination of the Cavity Fields 


By noting that the cavity field operator commutes with the atomic interaction Hamiltonian [A, iJint] = 0, then the 
Heisenberg equations of motion of the cavity field operators are determined by the single-particle Hamiltonian 

Eq. (A 6 ): dtAj = —i[Aj,HgQ^]/h— kAj, where the cavity-mode decay —k 4 i® included phenomenologically. They 
can be recast in the matrix form. 


4 (— I ( '^11 ~Q;i2\ / Ai\ ^ /?7i\ 

dt \A 2 / *^21 <322 J \A 2 / yl2) ' 


(Bl) 


where the elements of the ’’operator” matrix a are given by 


dll = (Aci -I- in) - 


A, 


2'^? ^ 

di 22 = (Ac2 + in) - 7 ^ All, 


di 2 — doi — 


(B2) 


If the cavity decay rate n is large, then the cavity fields reach steady states very quickly. By setting 9* Ai = 9 tA 2 = 0 
in Eq. (Bl), one can simultaneously solve the two equations of motion to obtain formal expressions for the steady-state 
field amplitudes Assj- However, one should take special care in solving these equations since the cavity fields and 
atomic fields commute with one another and this can give rise to ambiguities in solving these equations. In order to 
avoid such ambiguities, we symmetrize the equations of motion and exercise symmetrization procedure in all results 
following from the equations of motion. Thus, after setting dtAi = 9(^2 = 0 in Eq. (Bl), we re-express equations of 
motion as 


^diiAg 


-A. 


nan ) - 2 


^di2Ass2 


Ass 2 dl 2 


+ ^1 = 2 


^d22Ass2 + 4s2d22^ — ^ (q^2i4s1 + Assld21^ + ??2 — 0. 

(B3) 
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Equation (B3) can then be rearranged 

1 r 


^ssl — ^ + di2dii^) ^ss2 + ^ss2 (^ 116:12 + di2Q;i/) 

^ss2 — ^ (^22 ^21 “t“ d;2ld;22 ) ^Issl “1“ ^Igsl (^22 ^21 “t“ 0^210^22 ) 


+ idii^?7i, 


+ 'i'^22V2, 


(B4a) 

(B4b) 


where aii and q;^ 2 ^ are the inverse operators of ctu and q ;227 respectively, such that dii&ii = diidu = 1 and 
Q! 22 d^ 2 ^ = q;^ 2 ^q ;22 = 1 - In order to make the subsequent analyses somewhat easier and trackable, we assume that all 
dual variables (except rij at this moment) are equal, namely, Ai = A 2 = A, Ad = Ac 2 = Ac, and i^ae = = %■ 

We also introduce Ac = Ac + in for a shorthand. We expand the inverse operators to the second order in a small 
unitless parameter ^ = 2i^Q^/AAc ^ 1 (and with (Nr) ~ 10® one still has ^{Nr) ~ 10“^ ^ 1, see Sec. II for more 
details). 


^11 = I Ac - 


2 ^ 

A 


As 


-1 


-A 


-1 


« 2 - 2 ^ = ( Ac - ^ 


-1 


~ A: 


1 + 


1 + 


2 ^ 

AAc 

AAc 


A 2 + 


-Ai + 


4^ 

A 2 A 2 ' 

A 2 A? 


m 


N?,, 


(B5) 


such that diidd® = ai^dn = d 22 d^ 2 ^ = Q;^ 2 ^d 22 = 1 + Note that the error in symmetrizing Eq. (B4) is also 

of order 0(^®). This can be easily checked by substituting, say, Eq. (B4a) in the first equation of (B3). Equations 
(B4a) and (B4b) can now be simultaneously solved, yielding 


2 lssl — ff ^ 
2lss2 = if 


“t“ ^ (^11 0^120^22 CKr20^ri ^22 T 0^22 dr2 T ^22 ^120^11 ) 

;-i , ’ll fA-i, 


TI2Q.22 H-^ {^0.22 d.2lQ.ii + 0.210.22 + di^ 022 ^21 + d^i d2ld22) 


(B 6 ) 


where E = 


1 — (di 2 ci 2 i + d 2 idi 2 ) up to by noting q;i 2 = d^i = 2^g/A and (B5). We then have 


— 1 H-=— (di2d2i + ci2iQ!i2) — 1 H-=— (S-i-S— + , 

9 A 2 1 A2A2 V ^ 


2A? 


A 2 A? 


(B7) 


up to 0(C®). Using Eqs. (B2), (B5)-(B7), and retaining terms up to we obtain 


-^ssl — ~ 


^ss 2 — 


Ac 

i 

£ 


m 


2 %^ 

AAc 

2 ^ 

AAc 


(j]iN2 + II 2 S- 

(r]2Ni + 


+ 


4%" 

A 2 A? L 
4^ 
A 2 A? L 


V1N2 + f [S+S- + S-S+ ) + r^ 2 NS- 


^ +S.S+j+ m^S+ 


(B 8 ) 


By substituting steady-state solutions (B 8 ) and their Hermitian conjugates in the Hamiltonian Eq. (A 6 ), 

exercising symmetrization procedure again and retaining terms up to we can find an effective Hamiltonian which 
depends solely on the atomic operators. After some tedious though straightforward algebra, we obtain the cavity- 
field-eliminated effective many-body Hamiltonian 


iJeff = I d'^r [ ’5'^'Hso^ -I- ^ginj + ^52^2 + 9i2nin2 ) + 


)+ Y. + (U±££ + + 2UdsN£, (B9) 

^ r^l ,2 


where the cavity-field-eliminated, effective single-particle Hamiltonian density takes the familiar form 

"Hso = + kna^^ f] + Uext(r) + -k 

with effective two-photon detuning and Raman coupling given by 


A(A?+Ac 2)2 02 17lj, 


Hr = 


2 %^ 


A(A?-h«; 2)2 


A? - 

^ A 


H', 


1?1172 = 


R 


(A? -t k 2)2 


A? - 

A 


?7i?72- 


(BIO) 


(Bll) 
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(Note that 5' = 0, since we have assumed Aai = Aa 2 = Aa; cf. Eqs. (A2) and (A4).) The coefficients of the 
cavity-mediated long-range interactions are found to be 


^ A2(A2 + k2)3 "^2: 


_ 4fi%4A^(A2 - 3«2) ^ 

" A2(A2 + «;2)3 


4;i%4A (A2 _ 3^2) 

= A2(A2;.2)3 


U± = 


4h%4A, 


A2(A2 -hAc2)3 


- ivf + 


4h%^A 

^ A2(A2 


+ ^2)3 


- {V2 + 2?7?)k^ 


(B12) 


The terms with coefficients U 1 / 2 , and Uds in the effective Hamiltonian (B9) are the cavity-mediated long-range 

interactions. Note that in the special case of 771 = 772 = 77 , o ne has 6 = 0 and Ui = U 2 = U± = Uzp = Uds = U- _ 
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